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ABSTRACT 

We propose a simple estimator for the gravitational potential of cluster-size halos using the temperature 
and density profiles of the intracluster gas based on the assumptions of hydrostatic equilibrium and spherical 
symmetry. Using high resolution cosmological simulations of galaxy clusters, we show that the scaling 
relation between this estimator and the gravitational potential has a small intrinsic scatter of ~ 8% — 15%, 
and it is insensitive to baryon physics outside the cluster core. The slope and the normalization of the scaling 
relation vary weakly with redshift, and they are relatively independent of the choice of radial range used and 
the dynamical states of the clusters. The results presented here provide a possible way for using the cluster 
potential function as an alternative to the cluster mass function in constraining cosmology using galaxy clusters. 

Subject headings: cosmology: theory - galaxies: clusters: general - methods: numerical - X-rays: galaxies: 
clusters 



1. INTRODUCTION 

Clusters of galaxies are useful probes of cosmology. The 
evolution of their abundance probes the growth of structure 
and therefore p rovides constraints on c osmological param- 
eters (e.g., see IVoitl 12005b I Allen etal.ll2011l for reviews). 
Clusters are readily observable in multiple wavelengths. In 
addition to the optical, the regime where galaxy clusters 
are first observed, hot gas in the deep gravitational wells 
of clusters shines in X-ray. In the sub-millimeter regime, 
galaxy clusters appear as distortions in the cosmic microwave 
background (CMB) as the hot electrons in the intracluster 
medium (ICM) scatter off CMB photons through inverse 
Compt on scattering known as the Sunyaev-Zeldovich ( SZ) 
effect dSunvaev & Zeldovichlll970HCarlstrom et al.ll2002h . In 
recent years, cosmological constraints have be en yielded us- 
ing cluster counts alone in optical surveys (e. g. jBahcall et alJ 
20031 iGladders et alJl2007t iRozo et al.ll2010h. in X-ray (e.g., 
Vikhlinin et alJ l2009HMantz et alJ 120081 120101). and in SZ 
(e.g.. IVanderlinde et al] 120101: ISehgal et al] l20lil) . Cluster 
abundances and clustering also c an also give us co nstraints on 
primordial non-gauss ianity (e.g.JCunha et al. 2010) an d mo d- 
els of gravity (e.g.. lSchmidt et al.ll2009tlRapetti et alJl2010t) 

Ideally the simplest way to constrain cosmology is to mea- 
sure the abundance of observable signals from clusters and 
constrain cosmology directly from them, e.g., constraining 
cosmology through direct comparison of the cluster temper- 
ature function predicted from theory against observations. 
However, current theoretical models are unable to model clus- 
ters by their observables to the precision necessary for con- 
straining cosmological parameters down to percent level. On 
the other hand, theory does provide accur ate and precise pre- 
dictions of the cluster mass function (e .g., |jenkins et af]2001t 
I Warren ela l. 2006; Tinker et al. 20QH). The current paradigm 
of cluster cosmology therefore relies on estimating the clus- 
ter mass function from cluster observables via well-calibrated 
observable-mass relations. 

Instead of focusing on cluster mass, an alternative is to 
characterize clusters by their gravitational potential, a quan- 
tity which is also well predicted by theory. Gravitational 



potential has potentially three advantages over mass. First, 
observable quantities, such as gas temperature and gravita- 
tional lensing signals, depend on the gravitational potential 
directly and only indirectly on mass. Therefore gravitational 
potential provides a more direct connection to cluster observ- 
ables. Second, the shape of gravitational potential is more 
spherical than the matter distribution, and this may reduce the 
scatter in scaling relation due to the non-sphericity of mat- 
ter distribution. Third, different definitions of halo mass ex- 
ists, and they are not necessarily tied to the observable prop- 
erties of cl uster. For example, friends-of-friends (FoF) algo- 
rithm (e.g.. lEinasto et al.l [l984: Da vis et alJl 19851) . frequently 
used in defining halos in simulations, is difficult to imple- 
ment in observ ations. Also, FoF halo in general has irregu- 
lar shape (e.g., Lukic et al. 2009) and this poses another dif- 
ficulty when relating to observations, where in most cases 
measurements are made within some circular annulus. An- 
other way of defining a halo is through spherical overdensity, 
where average overdensity inside the halo is some value times 
the mean background density or the critical density of the uni- 
verse. While this definition is more observationally-oriented, 
the halo mass defined this way artificially evolves with red- 
shift as the mean background density and the critical density 
both change with redshift. 

Efforts have been recently made both theoretically and ob- 
servational ly in characterizing clus t ers w ith their gravitational 
potential. Angrick & Bartelmann (2009) derived the abun- 
dance of cluste r potential minima u sing a Gaussian random 
field formalism (Bar deen et al.ll 1986h . and was used to predict 
cluster temperature function using a simple spherical-collapse 
model to model the nonlinear e volution of the cluster poten- 
tial. On the observational side, IChurazov et all (120081 boiOh 
measured profiles of the gravitational potential of elliptical 
galaxies using both X-ray observations and stellar kinematics 
to constrain non-thermal pressure in those galaxies. It should 
then be straightforward, at least in principle, to extend these 
measurements to group and cluster-size systems. 

In this paper, we show that the cluster potential can be mea- 
sured reliably via observables, and that the cluster potential 
can serve as the defining quantity of cluster. We use Simula- 
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tions to estimate the gravitational potential of galaxy clusters 
by constructing a simple potential estimator based on intra- 
cluster gas prof i les, motivate d by the observational results of 
IChurazov etaTI (120081 120101) . In Section |2] we describe the 
simulations used in the paper. In Section [3] we describe our 
estimator of the cluster potential. In Section |4] we show the 
results of the scaling relations of this potential estimator and 
the gravitational potential. In Section[5]we give our summary 
and discussion. 

2. SIMULATIONS 

The simu l ation d ata we used are identical to those in 
Nagai et al. (2007a, b), where 16 cluster-sized systems are 
simulated using the Adaptive Refinement Tree jV-bo dy+gas- 
dynamics code ( lKravtsovlll999l:lKravtsov et~a l. 2002), which 
is an Eulerian code that uses adaptive refine ment in space 
and time, and non-adaptive refinement in mass (Kl vpin et ID 
1200 11) to achieve the dynamic ranges to resolve the cores 
of halos formed in self-consistent cosmological simulations. 
The simulations assume a flat ACDM model: f2 m = 1 - 51a = 
0.3, n b = 0.04286, h = 0.7 and cr 8 = 0.9, Q where the Hub- 
ble constant is defined as 100/i km s -1 Mpc -1 , and a% is the 
mass variance within spheres of radius 8/T 1 Mpc. The sim- 
ulations were run using a uniform 128 3 grid with eight lev- 
els of mesh refinement. The box size for CL101-CL107 is 
120 h~ l Mpc comoving on a side and is 80/T 1 Mpc comoving 
for CL3-CL24. This corresponds to peak spatial resolution of 
« lh~ l kpc and 5/T 1 kpc for the two box sizes respectively. 
Only the inner regions ~ 3-10/i _1 Mpc surrounding the clus- 
ter center were adaptively refined. The dark matter (DM) 
particle mass in the region around each cluster was m v ~ 
9.1 x 10 8 A _1 Af a for CL101-107 and m p ~ 2.7 x 1O 8 A _1 M 
for CL3-24, while other regions were simulated with a bot- 
tom mass resolution. 

In Table Q] we report rsooc ( tne radius of the cluster within 
which its average density is 500 times the critical density), 
Msoqc (the mass within rsooc), and our classification of the dy- 
namical state of the cluster based on the morphology of their 
m ock X-ray images. D etails of the classification can be found 
in iNagai et ail (l2007bl) . The cluster center is defined as the 
location of the potential minimum. 

We assess the effects of dissipative gas physics on the rela- 
tion between the cluster potential and gas properties by com- 
paring two sets of clusters simulated with the same initial con- 
ditions but with different prescription of gas physics. In the 
first set, the gas follows simple physics without any radia- 
tive cooling or star formation. We refer this set of clusters 
as non-radiative (NR) clusters. In the second set, metallicity- 
dependent radiative cooling, star formation, supernova feed- 
back and a uniform UV background were added. We refer this 
set of clusters as cooling+star formation (CSF) clusters. For 
a detailed description of the gas physics im plemented in the 
CSF clusters, please see Nag ai et all (12007a). 

3. ESTIMATOR OF THE GRAVITATIONAL POTENTIAL 

Following IChurazov et al.l d2008l) . we use a simple estima- 
tor of the cluster potential based on gas observables. Under 
the assumptions of hydrostatic equilibrium, the gradient of 

1 Note that this value of erg = 0.9 is higher than the current estimate of 
erg = 0.8. We do not expect the results and conclusion of the paper are going to 
change significantly. Nevertheless, a likely difference for a lower ag would be 
the later cluster formation time, leading to less dynamically relaxed clusters, 
potentially increasing the scatter of our potential estimator. 
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0.52 
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0.49 
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CL11 


0.90 


0.54 





CL14 


0.77 


0.51 
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0.35 
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the gravitational potential is related to the pressure gradient 
and density of gas by 



V0 = , 

Pg 



(1) 



where 4> is the gravitational potential, p g is the density of the 
intracluster gas, and P is the gas pressure. The gravitational 
potential can be obtained as a function of gas temperature T 
and gas density when we integrate Equation (Q]i by parts: 



dP 

Pg 



■ + constant; 



k B T 



k B T 
Limn 



d In p g + constant 
(2) 

where kg is the Boltzmann constant, /x = 0.59 is the mean 
molecular weight for the fully ionized ICM, and toh is the 
mass of the hydrogen atom. 

Physical properties relate not to the absolute value of the 
potential but the difference in the potential. The difference of 
the spherically averaged gravitational potential between two 
arbitrary radii r\ and r2 can be directly calculated from the 
temperature and density profiles. We define our potential dif- 
ference estimator which is based on the temperature and den- 
sity profiles of the intracluster gas as Aip, 



A^ = Tp(r 2 )-iP{r x ): 



kT(r) 



Limn 



kT(r) dlnpg 
/ioth dr 



dr. (3) 



One can note that iji is actually the enthalpy of the intracluster 
gas. If Aijj is a perfect estimator of the true potential dif- 
ference A<j>, then obviously we have A?/> = A<p, Since in real 
clusters, neither the gravitational potential is strictly spherical, 
nor the intracluster gas obeys perfect hydrostatic equilibrium, 
we expect Atp to deviate from A<f>. To assess how well Atp 
measures A<f> statistically, we fit a scaling relation of the two 
quantities using ordinary least squares: 



log Ac/) = a log (Atp/Aipo) +log A<j> , 



(4) 



where we set pmnAipQ = 20 keV. Assuming that the residu- 
als follow a log-normal distribution, we calculate the intrinsic 
scatter as the root mean square of the residuals <5i n A</> divided 
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by \/N—2 where N is the total number of clusters: 



2 

°lnAi/> ' 



N 



In A(j)i - a In ( Aip /Atpo) - In Acp Q ) 



(5) 

Note that we follow the convention of using natural log for 
the log-normal scatter. 

As Aip and Acp are dimensionally the same, we expect the 
slope of the scaling relation between them to be close to unity. 
We therefore also compute the scatter of the relation assuming 
unity slope, i.e., a = 1. In this case, we compute the intrinsic 
scatter as the root mean square of the residuals divided by 
y/N— 1 because of one less free parameter. We estimate the 
errors of the parameters a, log Acpo and a\ n Aip by generating 
10 5 bootstrap samples. 

We calculate the radial profiles for each cluster using log- 
arithmically spaced bins, with a total of 100 bins spanning a 
comoving radial range from 1 h~ l kpc to 5 h~ l Mpc. Note 
that our results are not affected by the exact details of the bin- 
ning scheme. The potential difference A(j> for each cluster is 
calculated from the potential profile directly measured from 
the simulation. The potential at each radial bin is calculated 
as the volume-average of the potential over all hydrodynamic 
cells inside the bin. 

To calculate Aip, we measure the spherically averaged 
gas density profile and temperature profile for each simu- 
lated cluster. The temperature profile is measured as mass- 
weighted temperature r mw which is the correct representation 
of the average specific gas thermal energy. For each radial 
shell, r mw is calculated as 



T = 

— -r-^ AT/' 



(6) 



where /?,-, 7), and AV, are the gas density, gas temperature, and 
volume of the cell i inside the radial shell and summed over 
all cells inside the radial bin. 

Although mass-weighted temperature can in princi- 
ple be me asured by combining X-ray and SZ ob serva- 
tions (e.g., lAmedio et all 120071 iMroczkowski et al.l 120071; 
lPuchwein"& Bartelmann 2007), X-ray observations alone do 
not measure the mass-weighted temperature. To approx- 
imate the temperature measured from X-ray spectra, we 
calculate Aip with t he spectroscopic-like temperature T s \ 
dMazzotta et al. 1 120041) which describes the temperature mea- 
sured by Chandra and XMM-Newton well: 



(7) 



This weighting scheme preferentially weighs more toward re- 
gions with high gas density and low temperature. In our calcu- 
lation of the spectroscopic-like temperature, we exclude cells 
that have temperature less than 10 6 K, which is well below 
the instrumental response of current X-ray instruments. We 
denote different Aip and the fitting parameters that use the 
mass-weighted temperature and the spectroscopic-like tem- 
perature by subscripts "mw" and "si", respectively. 

4. TESTING THE ESTIMATOR 

We use high-resolution cosmological cluster simulations to 
compare the potential difference estimator Aip with the real 
potential difference Acp directly measured from the simula- 
tion. We have made following assumptions for our estimator. 



Table 2 

Scaling Relation Parameters at z = 0.0 for ry = (0.0, 0.2,0. 4)r5oo c and 
n = r 500c 





Z = 0.0 


[l / r 500c , ^2 / ''500c 


] = 


[0.0,1.0] 


[0.2,1.0] 


[0.4,1.0] 




CSF 


0.896 ±0.098 


1.018 ±0.046 


1.056 ±0.066 




NR 


0.965 ±0.030 


0.970 ±0.035 


1.017 ±0.048 




CSF 


0.984 ±0.059 


1.010±0.043 


1.051 ±0.121 




NR 


0.951 ±0.041 


0.946 ± 0.054 


0.993 ±0.061 


log(/im H AflVmw) 


CSF 


1.607 ±0.022 


1.334±0.014 


1.362 ±0.033 




NR 


1.327 ±0.007 


1.333±0.014 


1.356 ±0.035 


log(/im H Acfo 1 si) 


CSF 


1.686±0.013 


1.375±0.015 


1.476 ±0.067 




NR 


1.336 ±0.008 


1.342 ±0.019 


1.370 ±0.037 


°ln A0,mw 


CSF 


0.1 82 ±0.027 


0.088 ±0.021 


0.140 ±0.035 




NR 


0.064 ±0.014 


0.085 ±0.019 


0.131 ±0.028 




CSF 


0.121 ±0.017 


0.072 ±0.016 


0.190 ±0.024 




NR 


0.075 ±0.010 


0.110±0.016 


0.147 ±0.026 




Fixed slope a = 1 




z = 0.0 


[nAsoociAsooc 


1 = 


[0.0,1.0] 


[0.2,1.0] 


[0.4,1.0] 


logt>m H A(/io ) mw) 


CSF 


1.597 ±0.020 


1.330 ±0.009 


1.335±0.015 




NR 


1.330 ±0.007 


1.341 ±0.009 


1.347 ±0.014 


log(/im H A0o,si) 


CSF 


1.686±0.013 


1.372 ±0.008 


1.419±0.019 




NR 


1.340 ±0.008 


1.360 ±0.012 


1.374±0.015 


°ln A0,mw 


CSF 


0.195 ±0.029 


0.091 ±0.023 


0.148 ±0.038 




NR 


0.070 ±0.015 


0.090 ±0.018 


0.1 36 ±0.029 




CSF 


0.125±0.017 


0.074 ±0.017 


0.198 ±0.020 




NR 


0.083 ±0.013 


0.119±0.0I9 


0.152 ±0.026 



1 . Input gas physics has little or no effect on the robustness 
of the estimator. 

2. The assumption of hydrostatic equilibrium is justified. 

3. The use of spectroscopic-like temperature is adequate 
to estimate the potential with little increase in the scat- 
ter of the A<p-Aip relation. 

4. The Acp-Aip relation does not evolve with redshift. 

5. The intrinsic scatter of the Acp-Aip relation follows 
log-normal distribution. 

In the following, we test the validity of the assumptions (1) 
- (4) and quantify the amount of scatter in the Acp-Aip rela- 
tion. For (5), we need a larger sample of clusters to quantify 
the deviation from log-normal behavior for the scatter which 
we leave it for future work. 

We report the values of the fitted slopes a, normalization 
loginmuAcpo) and intrinsic scatter <j\ n of the Acp-Aip scal- 
ing relations in upper halves of Table [2]- Table [4] In the lower 
halves of these tables, we report the normalization and intrin- 
sic scatter where we fit the Acp-Aip relation with the slope 
fixed to unity, a = 1 . We report values of the fitted parame- 
ters estimated using both mass-weighted and spectroscopic- 
like temperatures wherever possible. Since we are measuring 
and estimating the potential difference between the inner ra- 
dius r\ and the outer radius r^, for convenience we denote the 
radial interval for each scaling relation as [r\ , r{\- 

4.1. Effects of Dissipation 

Assuming hydrostatic equilibrium, we measure Acp and 
A^mw in the simulated clusters for both CSF and NR clus- 
ters at z = for different radial distance separations: [r\ , 7-2] 
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Figure 1. Top panels: scaling relations of the potential difference A0 and an estimator of the potential difference Aip constructed from gas density and mass- 
weighted temperature. The potential and the gas observables are measured from three sets of radial separations: from the cluster center to r^QQ c (left panel), from 
0.2r500r to r5oo c (middle panel), and from 0.4r5oo c to r^QQ c (right panel). The blue circles are for the CSF clusters and the red triangles are for the NR clusters. 
The blue solid line and the red dashed line are the fitted relation for the CSF and NR clusters respectively. The bottom panels show the corresponding residuals 
of the fits in the top panels. 



Table 3 

Scaling Relation Parameters at z = 0.0, 0.6. 1 .0 for r\ = 0.2r5oo<- and 



ln,r 2 ] = 


[0.2r 50Oc) r 5mc ] 


z = 


0.0 


0.6 


1.0 




CSF 


1.018±0.046 


0.936 ±0.050 


0.912 ±0.030 




NR 


0.970 ±0.035 


0.901 ±0.056 


0.916 ±0.054 


«sl 


CSF 


1.010 ±0.043 


0.852 ±0.070 


0.895 ±0.034 




NR 


0.946 ±0.054 


0.875 ±0.062 


0.881 ±0.069 


log(^m H A0 O ,raw) 


CSF 


1.334±0.014 


1.324±0.016 


1.314±0.018 




NR 


1.333±0.014 


1.300 ±0.025 


1.321 ±0.027 


log(£tm H A</>o iSl ) 


CSF 


1.375±0.015 


1.350 ±0.036 


1.380 ±0.023 




NR 


1.342±0.019 


1.313±0.030 


1.332 ±0.037 


ff ln A0,mw 


CSF 


0.088 ±0.021 


0.094 ±0.018 


0.069 ±0.011 




NR 


0.085 ±0.019 


0.150 ±0.022 


0.080 ±0.014 


""lnA^.sl 


CSF 


0.072 ±0.016 


0.150 ±0.022 


0.076 ±0.013 


NR 


0.110±0.016 


0.130±0.014 


0.107±0.014 




Fixed slope a = 1 


ln,n] = 


[0.2r 5mc .r 50 0c] 




0.0 


0.6 


1.0 


log(/im H A0 O ,mw) 


CSF 


1.330 ±0.009 


1.349 ±0.010 


1.353 ±0.009 




NR 


1.341 ±0.009 


1.345±0.013 


1.367 ±0.010 


log(/im H A</>o. s i) 


CSF 


1.372 ±0.008 


1.418±0.017 


1.436 ±0.010 


(a=l) 


NR 


1.360 ±0.012 


1.372±0.015 


1.401 ±0.014 


ff ln A0,mw 


CSF 


0.091 ±0.023 


0.102 ±0.021 


0.087 ±0.011 


NR 


0.090 ±0.018 


0.129 ±0.020 


0.099 ±0.026 


ClnA</>,sl 


CSF 


0.074 ±0.017 


0.172 ±0.024 


0.099 ±0.013 




NR 


0.119±0.019 


0.153 ±0.020 


0.136±0.035 



where we set r\ /rsooc = 0,0.2 and 0.4, and fix r2 = r^ooc, as it is 
the typical outermost radius current X-ray surveys are capable 
of measuring. This is also the ra dius within which the clus- 
ter is relaxed (lEvrard et al.lll996|) and whe r e the non-thermal 
pressure small, < 10% ("e.g.. lEvrardlll990t [Rasia et al] l2006t 
Nagai et al. 2007b). However, as we show later the estimator 
is robust to the exact choice of radial separation. 

In Figure [T] we plot the A0-A^ mw relations for 
[ri,r 2 ]/r 5 oo c = [0, 1], [0.2, 1], [0.4, 1] for both CSF and NR 



Table 4 

Scaling Relation Parameters at z = 0.0,0.6, 1.0 for r\ = 0.2r5oo c and 
r, = r 500c for A0 NR - Ai/> CSF Relation 



[n , r 2 ] = 


[0.2r5ooo>"500c] 


z = 


0.0 0.6 1.0 


«sl 


1.072 ±0.088 1.002 ± 0.074 1.024 ±0.068 
1.070± 0.084 0.867 ±0.122 0.997±0.086 


log(/tm H A0 o , raw ) 
log(^m H A0o. s i) 


1.271 ±0.024 1.288 ±0.022 1.278 ±0.036 
1.315 ±0.026 1.295 ±0.053 1.348 ±0.050 


°"ln A</>.mw 
CT lnA(/>,sl 


0.1 84 ±0.035 0.1 87 ±0.045 0.1 37 ±0.020 
0.164±0.031 0.270±0.054 0.162±0.030 




Fixed slope a = 1 


[n , r 2 ] = 


[0.2r5no c ,r 5 oo c ] 




0.0 0.6 1.0 


log(/tm H A0o,raw) 
log(^m H A0o. s i) 


1.255±0.019 1.287± 0.019 1.267 ±0.014 
1.297±0.017 1.356±0.028 1.350±0.017 


CHn A(/),mw 
^lnA^.sl 


0.1 82 ±0.036 0.1 80 ±0.042 0.1 33 ±0.029 
0.163±0.031 0.269±0.058 0.157±0.028 



clusters. All slopes of the relations are within the expected 
value of unity to within ler. The normalization log(/im H A7/io) 
is also near to the expected value of log20 « 1 .3. The intrin- 
sic scatter varies from ~ 6% to 15%, depending on the radial 
separation and input cluster physics. In general the resulted 
scatters of the fixed slope relations increase slightly compared 
to those where the slope is a free parameter. The normaliza- 
tions of the fixed slope relation remain unchanged to within 
lcr. The values of the fitted parameters are reported in Ta- 
ble 12 

For [r\ ,r2]/r5oor = [0.2, 1], [0.4, 1], the slope, normalization 
and scatter are similar for both CSF and NR clusters. Dissi- 
pative gas physics has little effect on the potential estimator 
outside the cluster core r\ > 0.2r^) C . 

For [ri , r2]/r5ooc = [0, 1], while the slope is similar in both 
CSF and NR clusters, the normalization of the CSF relation 
is offset to a higher v alue compared to the NR relation. As 
we discuss in Section l4~2l the higher normalization is due to 
strong gas rotation in the CSF clusters, where Ai/j underesti- 
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1 10 10 z 1 10 10 2 

Mm H A^ mw [keV] Atm H A^ mw [keV] 

Figure 2. Top panels: Scaling relations of the potential difference A<f> and 
an estimator of the potential difference Aip constructed from gas density and 
mass-weighted temperature. The potential and the gas observables are mea- 
sured from the cluster center to r^QQ c . The left panel shows the scaling relation 
without including pressure due to gas motions, while right panel shows the 
scaling relation with gas motions. The blue circles are for the CSF clusters 
and the red triangles are for the NR clusters. The blue solid line and the red 
dashed line are the fitted relation for the CSF and NR clusters respectively. 
The bottom panels show the corresponding residuals of the fits in the top 
panels. 

mates A</>, shifting the data points to the left in the A</>-A?/> mw 
plane. This underestimation is also shown in the upper left 
panel of Figure [3] where Aip mv/ /A(j) is plotted as a function 
of the inner radius r\ . The CSF relation have a relatively large 
scatter of about 18% due to strong dissipation in the cluster 
core regions. 

The potential estimator Aip is robust to dissipational gas 
physics outside the cluster core. As we show in Section [4721 
that once we correct for non-thermal pressure support due the 
gas motions, the estimator becomes robust to gas physics even 
when the cluster core is included. 

4.2. Non-thermal Pressure Support 

We test the validity of the assumption of hydrostatic equi- 
librium of our potential estimator by explicitly including the 
contribution of non-thermal pressure support by gas mo- 
tions, which are mainly responsible for the de parture of 
hydro static equilibrium in our simulated clusters (lLau et al.l 
2009). The effects of gas motions can be corrected by in- 
cluding the following term A-0 v>as in Equation (01, derived 
from the spherically sym metric radial Jean's equation (e.g., 
iBinnev & Tremainell2008l) : 

A^ = - (v*fo)-3(n)) -J^^dr-J" ^t^dr 

_ _ ' " (8) 

where and vj are the radial and tangential components of 
the mean-square gas velocities averaged over the spherical ra- 
dial shell, with respect to the cluster peculiar velocity defined 
as the average mass- weighted velocity of DM inside rsoo c . 

Figure|2]shows the scaling relation with and without the in- 
clusion of the non-thermal pressure term Atp v as . Including 
gas motions, i.e., replacing Aip mvl with Atjj mv/ + Aip v re- 
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Figure 3. Aip mvt /A<j> as a function of the inner radius r\ with the outer 
radius = rsooc at z = 0. The solid lines represent relaxed clusters and the 
dashed lines represent unrelaxed clusters. Top panels show the CSF clusters 
and the bottom panels show the NR clusters. The left panels show the fiducial 
m w without including gas motions, while the right panels include pressure 
support form gas motions. 

suits in better estimates of the potential difference by decreas- 
ing the scatter of the scaling relation to ~ 5% for [0, l]r5oo c 
for both NR and CSF clusters This is illustrated more clearly 
in Figure [3] where we plot the deviations of the estimated po- 
tential difference A-0 mw from the true potential difference A<fi 
as a function of r\ with r^ fixed at rsooe, with solid lines rep- 
resenting the dynamically relaxed clusters and dashed lines 
representing the unrelaxed clusters. Without the inclusion of 
gas motions, Atp mv , underestimates A<j> by about ~ 5%-20% 
for r\ > 0.1r5ooc in both NR and CSF clusters, a value that 
is consistent with gas motions biasing the ther mal pressure 
low by about the same amount (lLau et al.ll2009h . Most of the 
scatter is driven by the dynamically disturbed systems in the 
sample. For r\ approaching r^ = rsoo c , the scatter in both CSF 
and NR clusters increases as the outer regions of the cluster 
are less relaxed. In the inner region r\ < O.lrsooc, the devia- 
tions from hydrostatic equilibrium are approximately constant 
for the NR clusters, but for CSF clusters, A?/> mw increasingly 
underestimates Ac/) as r\ decreases due to strong gas rotation 
near the core w here the gas is rotationally supported induced 
by gas cooling dLau et al.ll201 ll) . Including gas motions takes 
into the account the effective potential due to the rotation, 
leading to a better recovery of the true potential difference 
and agreement between the scaling relations of the CSF and 
the NR clusters. This is shown in the upper-right panel of 
Figure |2] 

4.3. Spectroscopic-like Temperature 

The spectroscopic-like temperature T s \ (Equation (0) is 
generally less than the mass-weighted temperature T mw 
(Equation (O), as it weighs toward colder and denser gas. As 
shown in Table|2j for different [ri , rsooc] the scatter in the A0- 
A^i relation is also larger than A0-A-0 mw relation, as A-0 s i 
is biased towards the colder and denser gas, whose fractions 
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Figure 4. Top panels: scaling relations of the potential difference A</> and the potential difference estimator Ai/> mw for a fixed radial separation with the inner 
radius r\ = 0.2r5ooc and the outer radius ri = r^ c at different redshift: z = (left panel), z. = 0.6 (middle panel), and z = 1 (right panel). The blue circles are 
for the CSF clusters and the red triangles are for the NR clusters. The blue solid line and the red dashed line are the fitted relation for the CSF and NR clusters 
respectively. The bottom panels show the corresponding residuals of the fits in the top panels. 



and locations vary across different clusters. Nevertheless, the 
scatter is increased only by a few percent, with the exception 
for [ri,r2] = [0.4, \]r^) C where the large scatter is driven by 
a single cluster which has abnormally low T s \ due to a dense 
gas clump. 

Comparing the A0-A-0 s i relation between the CSF and NR 
clusters, we find that the normalization for the CSF relations 
are slightly higher than their NR counterparts at large r\ , be- 
cause of the relatively larger number of cold dense clumps in 
the CSF clusters than in the NR clusters. Nevertheless their 
normalization agree to within lcr. 

To account for the limited angular resolution of X-ray ob- 
servations, when we measure 7^, we average it over a radial 
window of Ar =100 kpc centered on the radius of interest. 
Varying the size of this window from 50 kpc to 200 kpc has 
little effect in changing the A<fi-Aip s \ relations. 

4.4. Evolution with Redshift 

Next, we investigate the evolution of the Aip-Aip rela- 
tions by fixing [n,r2] = [0.2, ljrsooc and fit the relations at 
z = 0.0,0.6 and 1.0. Figure [4] shows the resulting plots. All 
the relations at higher redshift are consistent with the z = 
results to within 2er. Scatter decreases slightly with increas- 
ing redshift and there is a weak trend of decreasing slope with 
redshift. The values of the parameters are shown in Table[3] 

The apparent weak evolution in the A<fi-Aip relation can 
perhaps be explained by the lack of change in the gravita- 
tional po tential over tim e once the DM halo has formed. For 
example. iLi et alj (120071) demonstrated semi-analytically that 
the circular velocity at the virial radius of the halo remains 
relatively unchanged throughout much of its mass accretion 
history. Given that our sample size is quite small (16 clus- 
ters), the apparent redshift evolution could be driven by a few 
clusters. A complete understanding of the redshift evolution 
of the potential well and a full characterization of the red- 
shift evolution of the Acp-Aip relation will require tests using 
simulations with more clusters and redshift outputs than the 
current study. 

4.5. Choice of Radial Separations 



In principle r\ and r^ that define A<j> = 4>(r2)-(f)(ri) and Alp 
can be chosen arbitrary. In practice, the best choice of [r\ , r2] 
would be the one that gives the lowest scatter. Once corrected 
for pressure from gas motions, the larger |n-r2| will result in 
less scatter, as the integral in At/j (Equation (O) encompasses 
larger radial range thus making the estimator more robust to 
small-scale variations in gas density and temperature. For ex- 
ample, setting [ri,r2] = [0.1,2]r5oo c gives a scatter of about 
~ 6%, compared to ~ 8% for [n,r2] = [0.1, ljrsooc- Decreas- 
ing r\ also gives smaller scatter as long as one stays away from 
the cluster core for the CSF clusters. For example,the scatter 
decreases from ~ 14% for [ri,r2] = [0.4, l]r5oo c to ~ 9% for 
[ r i) r 2] = [0-2, l]r5ooc- For r^ > rsooo the cluster gas deviates 
from hydrostatic equilibrium significantly with radius and Alp 
underestimates the true potential difference A(j>. However, 
this underestimation is small. As shown in Figure[5] the slope 
and normalization of the scaling relation are essentially the 
same as the radial interval changes from [n,r2] = [0.5, ljrsooc 
to [n,r2] = [0-l,2]r5ooc- The underestimation in the potential 
is significantly less than that of the hydrostatic mass estimate 
because At/j encompasses a large region of the cluster interior 
where deviation from hydrostatic equilibrium is small, while 
the hydrostatic cluster mass measured at a particular radius is 
dependent on the local value of the pressure gradient. 

Following similar logic, the choice of r\ and r^ should 
not be limited to functions of r^ c - Defining r\ and r^ 
through rsoor can be undesirable as r^o c depends on cosmol- 
ogy through the Hubble parameter H(z). One choice of the 
radial interval, independent of the halo radius, is simply the 
physical distance separations. Figure [6] shows the scaling re- 
lations for [ri,r 2 ] = [0.2, l]Mpc at z = 0.0, 0.6 and 1.0. For 
this radial separation the A<f>-Aip relation has a relatively 
small scatter of ~ 8%-9% for the given redshifts. The val- 
ues of slope and the normalization of the relations are similar 
to those with r\ and being functions of r^) C . The relation 
is the same for the different input gas physics. Redshift evo- 
lution of the relation is wea k but requires more investigation 
as discussed in Section l4"!4l 

4.6. Calibration for N-body Simulations 
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Figure 5. Top panels: scaling relations of the potential difference A</> and the potential difference estimator Ai/> mw at z = for inner radius ri = O.lrsooc and 
outer radius ri = 0.5f5ooc- (left panel), ri = rsoor (middle panel) and ri = 2r5oo c ( rl g nt panel). The blue circles are for the CSF clusters and the red triangles are 
for the NR clusters. The blue solid line and the red dashed line are the fitted relation for the CSF and NR clusters respectively. The bottom panels show the 
corresponding residuals of the fits in the top panels. 




Figure 6. Top panels: scaling relations of the potential difference A(f> and the potential difference estimator Ai/> mw for a fixed radial separation with the inner 
radius r\ = 0.2 Mpc and the outer radius ri = 1 Mpc at different redshift: z = (left panel), z = 0.6 (middle panel), and z = 1 (right panel). The blue circles are 
for the CSF clusters and the red triangles are for the NR clusters. The blue solid line and the red dashed line are the fitted relation for the CSF and NR clusters 
respectively. The bottom panels show the corresponding residuals of the fits in the top panels. 



Although becoming increasingly feasible, full cosmologi- 
cal hydrodynamical simulations of cluster formation with re- 
alistic physics like radiative cooling, star formation and active 
galactic nucleus feedback are expensive to run. A cheaper 
way would be using these simulations to calibrate less expen- 
sive dissipationless A^-body simulations. Since both pure N- 
body simulations and our NR clusters have gravity as the only 
driving physics, the results of our NR clusters should be sim- 
ilar to the that of the A^-body results. In this subsection, we 
therefore investigate how well our more realistic CSF clusters 
measure the gravitational potential wells of the NR clusters, 
and therefore the potential wells of clusters in A^-body sim- 
ulations. The results will be useful for calibrating a cluster 
potential function from pure A^-body simulations. 

We estimate the NR potential difference A(/> NR using the 



potential estimator of the CSF clusters Ai/>^ F . In Figure [JJ 
we show the comparison between the A0 NR -A^^ F relation 
for [ri,r2] = [0.2, ljrsooc at z = 0.0,0.6, and 1.0 (in blue solid 
lines) and the A0 csf -At/>^ f relation (in black dashed lines). 
The slope of the A</> nr -At/>^ f relation is steeper than the 
A0 csf -At/i^ f relation. This is because dissipation results in 
deeper potential wells in the CSF clusters compared to their 
NR counterparts, and this effect is stronger in small size halos 
which have lower virial temperature. Therefore Ai/jJ^f over- 
estimates the potential difference A0 NR more for the smaller 
size systems, leading to a steeper slope in the A</> NR - Aip^ 
relation. The scatter in A0 NR -A-0^ F relation also increases 
because the NR clusters are not in exact identical dynamical 
states as their CSF counterparts as the dynamical evolution 
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Figure 7. Top panels: scaling relations of the potential difference A</> NR in the NR clusters and the potential difference estimator Ai/jJJj^f of the corresponding 
CSF clusters, for a fixed radial separation with the inner radius r\ = 0.2r5oo c and the outer radius r% = 0.5r5ooc at different redshift: z = (left panel), z = 0.6 
(middle panel) and z = 1 (right panel). The blue circles are for the CSF clusters and the red triangles are for the NR clusters. The blue solid line is fitted relation. 
The black dashed line is the best fit for the A0 CSF -A?/£ F relation. The bottom panels show the corresponding residuals of the fits in the top panels. 



of each cluster halo is affected by the dissipative gas physics 
which changes the structure and hence dynamics of the halos. 
Table |4] summarizes the fitted parameters. 

5. SUMMARY AND DISCUSSION 
5.1. Summary of Key Results 

In this paper we propose a simple estimator for the gravi- 
tational potential difference in clusters of galaxies. This esti- 
mator is based on the density and temperature profiles of the 
intracluster gas under the assumptions of hydrostatic equilib- 
rium and spherical symmetry. 

Using high resolution cosmological hydrodynamical simu- 
lations of galaxy clusters, we have tested that this estimator 
is a robust and accurate estimator of the true gravitational po- 
tential of simulated clusters. The key results are summarized 
below: 

1. The scaling relation between this estimator and the true 
gravitational potential difference A(j)-Atfj has an intrin- 
sic scatter of — 8%-15%. 

2. Input gas physics has little effect on the Atp-Atfj rela- 
tion as long as the cluster core is excised. 

3. The assumption of hydrostatic equilibrium is justified 
outside the cluster core. Adding non-thermal pressure 
support reduces the scatter in the A<j>-Aip relation for 
all radii. The core need not be excised when non- 
thermal pressure support is included. 

4. With the core excised or kinetic pressure from gas mo- 
tions included, the slope and normalization of the Acj>- 
Aijj scaling relation are relatively independent of the 
choice of radial separations. The scatter of the relation 
decreases as the separation increases. 

5. The use of spectroscopic-like temperature is adequate 
to estimate the potential with little increase in scatter in 
the relation. 

6. The A(fi-Aip relation evolves weakly with redshift. A 
full characterization of redshift evolution is needed. 



The results presented above suggest that gravitational po- 
tential can serve as a useful alternative to cluster mass for 
cluster cosmology. X-ray mass-observable scaling relations 
are usually dependent on baryonic physics. For example, the 
normalization of the current tightest mass- observable relation, 
the M-Yx relation (Kra vtsov et alj 120061) . is dependent on 
non-thermal pressure support; the normalization of the cluster 
mass-ga s mass M-M pas re lation, is affected by the input gas 
physics dFabjan et al.ll201 II) . Our results show that the poten- 
tial scaling relation, with appropriate radial range selection, is 
little affected by both of these effects. 

Gravitational potential has another advantage over mass in 
that it is more spherical than that of the matter distributionQ 
The assumption of spherical symmetry works better for grav- 
itational potential than matter density, and this is perhaps 
partly why the potential scaling relation has relatively small 
scatter compared to mass-observable relations. 

An essential characteristic of the gravitational potential 
scaling relation is that we are free to choose the radius where 
we measure the potential, not limiting to the halo radius de- 
fined in terms of overdensity with respect to the mean den- 
sity of critical density of the universe, both of them functions 
of cosmology. This takes away the rather arbitrary nature of 
cluster mass, where different cluste r mass definitions give dif- 
feren t cluster mass functions (e.g., White 2002; Tinke r et"aT] 
2008). 

Despite the advantages given above, there are also several 
caveats and uncertainties in using the proposed potential esti- 
mator for cluster cosmology. We see a weakly evolving red- 
shift evolution in the potential scaling relation that needs to be 
understood well before using gravitational potential for clus- 
ter cosmology. Although the potential measured at different 
radial separations result in essentially the same potential scal- 
ing relation, it is uncertain whether the different radial sepa- 
ration would result in different cluster potential function. Fur- 
thermore, we have assumed spherical symmetry for the grav- 

2 We give a comparison of the shape of the dark matter distribution and 
the shape of th e gravitational p otential well of the same set of clusters used 
in this paper in Lau et al. 1 201 1). 
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itational potential and the gas distribution, but have not tested 
for the effect of deviation from spherical symmetry, although 
we expect that the effect to be small compared to the effect on 
cluster mass. Since our gravitational potential estimator in- 
volves an integral over radius of gas temperature and density, 
which is the dominant term in the estimator (Equation (01), 
the gas density and temperature profiles may need to be mea- 
sured with high spatial resolution. Accurate measurements of 
gas temperature also require high-resolution X-ray spectrom- 
eters. This can be technically challenging, especially if we 
want to measure the potential of high-z clusters. However, we 
note that combined X-ray/SZ analysis is able to recover the 
temperature prof ile without the need of X-ray spectra (e.g., 
iNord et alJl2009h . 

5.2. Toward the Cluster Potential Function 

Given the results presented in the paper, it will be interest- 
ing to see whether the cluster potential function is indeed a 
better alternative for constraining cosmology than the com- 
monly adopted cluster mass function. In this subsection, we 
discuss ways to improve the potential scaling relations and to 
construct the cluster potential function. 

We suggest that future modeling efforts toward the goal 
of using the cluster potential function for cluster cosmology 
should focus on two main areas: (1) further investigation on 
the intrinsic nature of the potential scaling relation and the 
properties of the cluster potential function and (2) implemen- 
tation of the potential scaling relation in actual observations. 

The first aspect aims to answer the question whether the 
cluster potential function can replace the cluster mass func- 
tion for cosmology. To see whether the cluster potential func- 
tion has more constraining power than the cluster mass func- 
tion, comparison between the cluster potential function and 
the cluster mass function calibrated from the same large-scale 
simulation is required. The comparison will have to focus on 
whether cluster function recovers the fiducial cosmological 
parameters better in terms of degeneracies between param- 
eters, statistical errors, and systematic uncertainties. Simu- 
lations with statistically large sample of clusters at multiple 
redshifts will help constraining the weak redshift evolution of 
the potential scaling relation shown in the current work. They 
will also help further characterize the scatter and systematics 
of the potential scaling relation, for example, deviations from 
lognormal distribution in the scatter. 

The second aspect focuses on developing new ways of 
characterizing gravitational potentials based on cluster ob- 
servables and apply them to synthetic and real observations. 
For example, the potential estimator proposed in the cur- 
rent work will benefit from mock X-ra y analyses in address- 
ing systematics like gas clumping dMathiesen et al.l 11999k 
Simio nescu et aljfeoi It iNagai &Laull201lh and some of its 
observational challenges discussed in the previous subsec- 
tion. It is reasonable to expect that other types of cluster 
potential estimators exist and may even perform better than 
the one proposed in this paper. For example, gravitational 
lensing signals like shear and convergence depend on the pro- 
jected gravitational potential and may provide a new probe 
to the gravitational potential. Velocity dispersion of cluster 
galaxies directly probes the cluster potential, although it can 
be affected by velocity anisotropy and projection effects. SZ 
signals measure the integrated gas pressure, and it is similar 
to the gas enthalpy-based estimator proposed in the current 
work. Multiple probes of the gravitational field will be help- 
ful in determining systematics, and the covariance of differ- 



ent potential estimator may also improve the precision of the 
measured potential , as it is for the case of cluster mass proxies 
dStanek et alj2010l) . 

Clearly, much more work is needed to address these issues 
and we hope that this paper provide a starting point for the 
investigation of using the cluster potential function for cos- 
mology. 
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